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We study the dynamical evolution of perturbations in the gravitational field of a collapsing fiuid 
star. Specifically, we consider the initial value problem for a massless scalar field in a spacetime 
similar to the Oppenheimer- Snyder collapse model, and numerically evolve in time the relevant wave 
equation. Our main objective is to examine whether the phenomenon of parametric amplification, 
known to be responsible for the strong amplification of primordial perturbations in the expanding 
Universe, can efficiently operate during gravitational collapse. Although the time-varying gravita- 
tional field inside the star can, in principle, support such a process, we nevertheless find that the 
perturbing field escapes from the star too early for amplification to become significant. To put 
an upper limit in the efficiency of the amplification mechanism (for a scalar field) we furthermore 
consider the case of perturbations trapped inside the star for the entire duration of the collapse. In 
this extreme case, the field energy is typically amplified at the level ~ 1% when the star is about to 
' cross its Schwarszchild radius. Significant amplification is observed at later stages when the star has 

even smaller radius. Therefore, the conclusion emerging from our simple model is that parametric 
amplification is unlikely to be of significance during gravitational collapse. Further work, based on 
CNJ ' more realistic collapse models, is required in order to fully assess the astrophysical importance of 

\ parametric amplification. 
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I. INTRODUCTION 
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. A. A brief background 
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The gravitational collapse of stellar objects with mass over a certain threshold is known to produce one of the 
most spectacular events in Universe. Due to the release of an enormous amount of energy, collapse is one of most 
favourable sources for all kinds of electromagnetic observations. It is expected (or hoped) that the same will come 
true for gravitational wave observations, which are about to begin in the next few years with the completion of a 



I network of interferometric detectors (for a recent review see [Q). A considerable amount of theoretical work (numerical 



in its majority) has been devoted to the description/prediction of the emitted gravitational waveform produced by 
• • , such a catastrophic event (see for instance ||]). Despite this effort, great uncertainties still remain, mainly due to 
the unknown initial state of the collapsing body, the complexity of the physics involved in realistic collapse and the 
related excessive computational requirements. 
• The simplest general relativistic collapse model was introduced more that sixty years ago with the seminal work 
■ of Oppenheimer and Snyder (hereafter 0-S) ||^. In that model the collapsing star is approximated as a freely falling 
"dust" fluid ball (see Section II for more details). Using the 0-S model. Price Q was the first to study the evolution 
of scalar field perturbations during spherically symmetric gravitational collapse, focusing on the last stages where 
a black hole is formed. Subsequent studies by Cunningham, Price and Moncrief |^ were extended to gravitational 
perturbations. It was found that, at the final stage of the collapse, the emitted waveform was dominated by the 
exponentially decaying quasinormal mode "ringing" of the newly born black hole. At even later times, the signal 
decayed as a power-law, representing radiation backscattered from the asymptotic gravitational field. 

Collapse of relativistic polytropes was modelled by Seidel Q (for a non-rotating configuration) and by Stark and 
Piran (for an axisymmetric rotating configuration) , who provided an estimate of the amount of energy released via 
gravitational radiation. In their calculations, the resulting waveform was mainly comprised of a short burst followed 
by the black hole ringdown. Recent efforts include the relativistic rotational core-collapse study of Dimmelmeier et.al. 
1^ and the more general investigation by Fryer et.al. ||]. We also refer the reader to the article of Font |l^] (and 
references therein) where numerical techniques in general relativistic hydrodynamics are discussed. 
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B. Motivation 



In the present paper we aim to study a specific aspect of gravitational collapse that (to the best of our knowledge) 
has been overlooked in the literature, namely, the possibility of having parametrically amplified perturbations within 
the star. That such a process could operate, at least in principle, can be justified with the following argument. 

It is well known that there is a qualitive similarity between the time-dependent gravitational field of a spherically 
symmetric collapsing star and the gravitational field of the closed Friedmann cosmological model. This similarity has 
been quantified in the celebrated 0-S collapse model which, despite its simplicity, gives reliable results. On the 
other hand, one of the most important results in the study of cosmological perturbations is that their interaction with 
the background gravitational field of the expanding Universe can lead to considerable amplification (in addition 
to the overall "adiabatic" change ^ a~^, where a is the cosmological scale factor). This mechanism is believed to 
be responsible for the immense amplification of tiny primordial gravitational quantum vacuum fluctuations and their 
pro motion to classical fields, which are likely to have left measurable imprints in the CMB radiation observed today 
P^ . Moreover, these fields put a strong candidature for direct observation by the gravitational wave detectors 

Hence, based on the two points above, a natural question that arises is whether the gravitational field of a collapsing 
star could pump energy into perturbations (that can be of fluid or gravitational nature) temporary "living" inside 
the star. In this work we plan to give an answer to this question by considering the simplest possible model for a 
collapsing star. In particular, we keep those pieces of physics which are related to the phenomenon we try to study, 
and strip off unrelated features which might be a cause of technical problems and confusion. We shall adopt a collapse 
model which very closely resembles the known 0-S model, but is even more simplified regarding the exterior to the 
star gravitational field (see next Section). In a realistic situation, we can think of a collapsing star as a generator of all 
kinds of perturbations (fluid, gravitational). In our toy model, this process is "mimiced" by placing some initial field 
inside the star and studyingits subsequent evolution. Moreover, we are free to choose the time at which this field is first 
placed inside the star. Finally, another major simplification is to consider a massless scalar field as the perturbation. 
We expect, however, that our results can be extrapolated to the realistic case of gravitational perturbations. This 
claim is based on the work of Ford and Parker jT^ , where both scalar and gravitational perturbations in Friedmann 
spacetime were considered. 

The remainder of the paper is organised as follows. In Section II we give a brief review of the 0-S collapse model 
(subsection IIA) and subsequently discuss the properties of the scalar wave equation in the field of a collapsing star and 
of full Friedmann spacetime (subsection IIB). In subsection IIC we introduce our simplified collapse model. Section III 
is devoted to our numerical results (time-evolutions). Analytical calculations regarding the amount of amplification 
in terms of the field's energy, are presented in Section IV. Section V offers a physical insight into the results of the 
previous Sections. Finally, a concluding discussion can be found in Section VI. Two Appendices are devoted to some 
technical details. Throughout the paper we have adopted geometrised units c = G = 1. 



II. A TOY-MODEL FOR GRAVITATIONAL COLLAPSE 



A. The Oppenheimer-Snyder model 

We set off by giving a brief review of the standard 0-S model (for a detailed presentation see |l^]). In this model, 
the spacetime of a collapsing, spherically symmetric, pressureless, homogeneous fiuid "ball" is constructed by patching 
together a piece of Schwarzschild geometry (describing the vacuum outside the star) and a piece of closed (fc ~ +1) 
Friedmann geometry (describing the stellar interior). The two metrics are matched smoothly across the surface of the 
star. Explicitly, the interior metric is written in familiar comoving coordinates (?], x, ^, 0), 

ds^ = a^{r])[~dr]^ + dx" + sin^ x{dO^ + sin^ Odcf^)] (1) 

while the exterior metric is, in Schwarzschild coordinates, 

ds^ = ~f{r)dt^ + f-\r)dr^ + r^{d9^ + sin^ 0d,p^) , (2) 

where /(r) = 1 — 2M/r. The scale factor in the Friedmann domain is, 

a{v) = ^a„,{l + cosjj) (3) 

with Om denoting its initial value (which is also the maximum value). The surface of the star is labelled as Xo and R in 
these two coordinate frames and as seen in Schwarzschild coordinates is moving inwards as a timelike radial geodesic 
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with zero initial velocity at some distance i?i corresponding to the initial stellar radius at the onset of collapse. This 
translates into 

i?(j?) = iiii(l + cosr?) . (4) 
The relation between the two different time coordinates is given by 



t = 2M\n 



(i?i/2M- 1)1/2 +tan(r//2) 



(iii/2M- 1)1/2 _tan(77/2) 



+ 2M{Ri - 1)1/2 ^ {Ri/AM){ri + sin??)] . (5) 



Note that the conformal time coordinate rj is restricted between r] = (initial time when collapse begins) and rj = w 
at which moment the star has shrank to zero radius. Of course, the physically relevant value rjc for r] is the time when 
the star is crossing its Schwarzschild radius, 

rjc = arccos — 1 . (6) 



The initial stellar radius as measured in units of the stellar total mass M is, 

^ = 2sin-2xo. (7) 

The parameter Xo can take any value between {Ri/M — » oo) and 7r/2 {Ri/M = 2). A typical collapsing object has 
Ri ~ 10^ — lO^M (the size of a white dwarf), and this value can rise even higher for supermassive stars. 



B. The scalar wave equation in the field of a collapsing star 

As we mentioned earlier in this paper, we shall consider the simplest possible case of perturbations, i.e. a massless 
scalar field which obeys the covariant wave equation V"Vq;$ = 0. We can immediately separate the angular 
dependence by means of the decomposition 

^ = J2RemYim{0,cl>) ■ (8) 
In the Friedmann domain (also simply termed "interior" hereafter), the wavefunction R^m can be rescaled as 

O / X '0(?m(??,X) /OX 

a{V) sinx 

Note that by this operation wc have isolated the usual "adiabatic" change ^ a~i due to the contracting background. 
The new wavefunction tp satisfies (for brevity we drop the £, m subscript hereafter), 

dli,-dl^ + v{r^,x)^ = Q , (10) 

where the effective potential V is given by (for x < Xo), 

V{71,X) = . 2 - 1 - - = ■ 2 - COS 2(^/2) . (11) 

sm X a sin X z 

Similarly, for the Schwarzschild domain (to be called "exterior") we have, 

R{t,r) = ^ip{t,r) (12) 

with tj) satisiying the familiar Regge- Wheeler equation, 

dl^l) - d^^ij - F(r)V = (13) 

with 
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(14) 



and where r* is the so-called "tortoise" radial coordinate defined by dr^:/dr = 1/ f{r). 

Naturally, our attention will be mainly focused on the interior time-dependent field which is identical (by construc- 
tion) to the field of a closed Friedmann model. The exterior field is of little concern to us, as it is responsible for 
effects like the quasinormal mode ringing and late time power-law tails which are the expected (and well studied) 
dominant components of the emitted signal at the final stages of a realistic collapse scenario. 

As a warm up, we first consider the initial value problem for the wave equation ( p^ ) in a pure Friedmann gravitational 
field. Moreover, we restrict ourselves to the monopole case £ = 0. Assuming static initial data ih(Q,x) = '0o(x)j 
d,fip{0, x) — ^'^d imposing the additional boundary condition '0('7j '"') = 0, the general solution of (nOh is given by. 



+ 00 r 



V'('?,X) = -y]sin(nx) / 
^t^, Jo 



dx V'o(x) sin(nx) 



cosinrj) H ta,n(r]/2) sininrj) 

2n 



(15) 



The key feature of eqn. (^_5|) is the term ~ tan(?7/2). This term originates from the d/a term in the potential eqn. 
( prj ) and consequently encodes all the effects associated with the time variance of the background gravitational field. 
We can apply the solution (^5|) for initial data of the form of a narrow Gaussian pulse. The resulting field is shown 
in Fig. |l|. 
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FIG. 1. Scalar field propagation in Friedmann spacetime. The initial field is shown on the left and later "snapshots" are 
shown on the center of the figure {rj — OAn) and on the right {rj = O.SStt). 



The field initially evolves as in flat spacetime, but at later times, we clearly observe a growth in the amplitude 
which is more pronounced in the "wake" behind the travelling pulse. This wake is mainly comprised of low-frequency 
waves (which "feel" more strongly the background gravitational held) as opposed to the main pulse which includes 
the higher frequency contribution. As 77 approaches tt the field becomes enormous as compared to its initial amplitude 
and eventually diverges, as predicted by eqn. (p^. This behaviour is a typical example of parametric amplification 
in Friedmann spacetime. 



C. Pushing the simpHfication further 



Let us now return to the 0-S model. The difference in the field evolution from the preceding picture is obvious: 
the field, as it propagates outwards, will reach the stellar surface at Xo- Once there, it will partially transmit to 
the exterior and partially reflect back to the interior (see Fig. |^). In principle, one could try and describe such a 
process analytically, by making some simplifications for the exterior field (for example, by keeping only the centrifugal 
term in the potential). However, such an attempt soon runs into serious technical problems related to the coordinate 
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"discontinuity" at the stellar surface. To overcome this difficulty, one could use comoving Novikov coordinates 
which smoothly covers the entire spacetime. However, in this case the resulting exterior wave equation is not separable, 
an undesired feature for any analytical work. At the same time, we should be almost certain that this coordinate- 
induced problems has nothing to do with the physics we try to study here. 



exterior 




XO X 

FIG. 2. A schematic illustration of our problem: Initial data placed inside a "collapsing" star propagate inwards and outwards 
experiencing partial and total reflections at the stellar surface (x ~ Xo) and at the center (x ~ 0) respectively. 

Of particular relevance to this point is the recent study by Laarakkers and Poisson fl^ ] on scalar wave propagation 
in the so-called Einstein-Strauss spacetime. This spacetime is, as the 0-S model, a combination of Schwarzschild 
and Friedmann geometries but patched in the opposite sense with the interior being Schwarzschild and the exterior 
Friedmann. Apart from providing fully numerical results, these authors considered a simple toy-equation that clo sely 
resembles the pair of exact wave equations (p^, jl^). The construction was based on simply replacing eqns. (p^,(|l^) 
with a single wave equation by extending the Friedmann coordinates over the entire spacetime. It was shown that the 
simplified wave equation captures most of the effects present in wave propagation in the Einstein-Strauss spacetime. 

For the purposes of the present paper we have chosen to follow a similar path and replace eqns. (p^),(p^ with the 
toy-equation, 



where the new effective potential is, 



y(r/,x)V' = 



(16) 



for X > Xo, exterior 
icos~^(r7/2) for X < Xo, interior 



(17) 



Hereafter, we shall consider only the monopole case £ = 0. We have verified that all results presented in this paper 
remain qualitively the same even for higher multipoles. The important point is that the new potential is identical 
to the "exact" 0-S potential in the interior, where the phenomenon of parametric amplification would take place. 
On the other hand, our model is unable to describe any effects related to the exterior gravitational field, such as 
backscattering of outgoing waves. In principle, this last feature could be of importance, as it would lead to trapping 
of some portion of the initial field inside the star. For this reason, we have also considered a star whose surface fully 
reflects any outgoing waves (see Section IV). 

In a sense, eqn. ( p^ describes a static star (as seen by both external and internal observers) with some kind of 
time-varying gravitational potential in its interior. This star, and with respect to external observers only, knows 
nothing about collapse. Physical parameters like the total duration of collapse have to be "borrowed" from the 0-S 
model. We believe that the wave equation (|T^), just as in the case of Einstein-Strauss spacetime, yields reliable 
predictions on scalar wave dynamics. 

It is possible to find an analytic solution to the above toy-equation (see [|l^ for a treatment of the fc = Friedmann 
equivalent) that describes the first transmission/reflection of a wavepacket through the potential discontinuity. In 
principle, we could extend this solution to incorporate multiple transmissions/reflections for the needs of our problem. 
However, the resulting expression becomes increasingly messy. This is the main reason that led us to approach the 
problem numerically. 
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III. TIME EVOLUTIONS 

We have written a simple numerical code for solving eqn. ( p^ ) after prescribing some initial data (typically in the 
form of a Gaussian pulse). The code, which employs a standard leapfrog step algorithm, was tested for convergence 
and stability. We have widely experimented with collapsing models of different i?i and placing the initial data at 
different stages of the collapse. A typical example of our results for the field as observed outside the star is shown 
in Fig. ^. In this simulation we placed the initial data at = 0, i.e. at the onset of the collapse, and centered at 
X = Xo/2, while the initial radius is Ri = W^M. We see that the field has escaped from the star, essentially unaffected, 
after a lapse of time of the order ~ 3xo/2 (which is the light time-travel for the, initially ingoing, wave component 
to reach the surface) . Only a minute fraction of the initial field (specifically the low frequency components) has been 
left inside the star while R ^ M. The same scenario is repeated in Fig. ^, but now the initial data is placed at a 
later time rji — 0.977r. This corresponds to an instantaneous stellar radius R « 22M. The final time is rjf « 0.997r 
which gives a radius R w 2.1M. In this case, the potential inside the star is much more pronounced and consequently 
is more effective in trapping inside the star a small fraction of the initial pulse for a considerable time. The majority 
of the field, however, escapes during the first two crossings of the stellar surface. In the same Figure we illustrate the 
signal from a star with a potential "frozen" at its value aX ij = rji (dashed curve) and at 77 = 77/ . We can see that the 
general appearance of the signal is similar for a static and non-static potential. It is very hard to say if there are any 
genuine features of parametric amplification present in the signal. For this reason we illustrate (see Fig. ^) the field 
as it looks at 77 = 0.9987r when the stellar radius has decreased to i? « O.OOIM. An amplified field inside the star 
can be clearly seen now (although it could never escape to infinity). In all the situations we examined, we reached a 
similar conclusion. Amplification becomes notable only after the star has crossed the Schwarzschild radius. Our final 
time evolution concerns the appearance of the emitted field as a function of the frequency content of the initial data. 
This can be simply done by considering as initial data a modulated Gaussian pulse, 

Mx) = exp[-0.5(x - Xc)V^'] sin[A(x - Xc)] , (18) 

where Xc the pulse location and A is the modulation frequency. In Fig. |^ we show a comparison between the signals 
originated by an unmodulated pulse and a pulse with A ~ 20/xo- AH other parameters are the same to the ones of 
Fig. ^. As we should have expected, the high-frequency pulse is escaping from the star more easily, and the late time 
"tail" is strongly suppressed. 
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FIG. 3. The field as observed at the radial location x = l-5xo, for Ri — W^M. The initial pulse was placed at Xo/2 and 
r? = 0. 
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FIG. 4. The field (soUd curve) as observed at x = l-5Xo, for Ri = lO^M. The initial pulse was placed at rji = 0.977r 
(corresponding to a radius R « 22M). The dashed and dotted curves depict the field for a "frozen" static potential V{ri = r]i) 
and V{r] = rjf) respectively. 
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FIG. 5. "Snapshots" of the field at time rj = 0.99 987r (solid curve) corresponding to a radius R ~ O.OOIM and at r; = 0.997r 
(dashed curve) corresponding to R~ 2M. In the former, it is easy to distinguish the amplified field in the stellar interior. The 
initial radius is Ri = 10*M, and the initial data were placed at rn = 0.987r. 
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FIG. 6. Dependence of the observed signal (at x = 1.5xo) on the frequency content of the initial data: the dashed (solid) 
curve was generated by a "simple" (modulated) Gaussian pulse as initial data. It is obvious that high frequency components 
escape from the star much easier. 



IV. A "MIRROR" STAR 



The results of the previous Section are clearly not very optimistic regarding the efficiency of parametric amphfication. 
In order to put an upper limit to this efficiency (at least for scalar perturbations) we have considered a "mirror" star, 
i.e. a star which totally reflects any outgoing field at the surface. Clearly, this situation is the most favoured one 
from the point of view of the amplification process, as the field will stay trapped inside the star for the entire collapse 
event. 

Mathematically, this mirror star is modelled by setting the additional boundary condition V''('7j Xo) — 0. Hence, the 
field inside the star has two inflection points at x = 0: Xo- It is straightforward to find an analytic solution to this 
problem that also satisfies static initial data V'(0, x) = V'o(x) (see Appendix A for details), 

2 

V'('7, X) = — y] C'n sin(n7rx/Xo) 

where 

Cn= f dxipo{x)sin{mrx/Xo) ■ (20) 
Jo 

We can strictly quantify the amplification of the field by means of the total energy 

E=ll^" dx[{d,ijr + {d^^f + V^'] . (21) 

It can be easily shown (see Appendix B) that 

E{ij) - E{0) = \ f r dij'dx ^^d,yV . (22) 
^ Jo Jo 

In this simple way we can see that the field energy is not conserved due to the time-varying potential. 

Combining eqn. (|l9|) with eqn. (p2[), we can monitor the energy of the scalar field inside the collapsing star. Results 
of this calculation, for a selection of i?;, are shown in Fig. ^ where we plot the quantity logj^Q |100x {£^(77) — i?(0)}/i?(0)| 
(i.e. the percentage fractional change in the energy) as a function of 77. The conclusion from this calculation is rather 
dissapointing from an "observational" point of view: The amplification factor is ^ 1% when the star is about to 



cos(7i7r?7/xo) 



Xo 

- — tan(7;/2) sm{mn]/xo) 



(19) 
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cross its Schwarzschild radius. It only becomes important, at later times, when the radius is close to zero. This is in 
agreement with the behaviour seen in the time-evolutions in the field of a "transparent" star (Section III). According 
to Fig. ^ there is a ~ 10% amplification when R 0.25M. It is interesting to note that, according to Fig. ^ the 
amplification factor is almost the same for a given value of R and different Ri. An intuitive explanation of this 
behaviour is given in the following Section. It is important to mention that the amplification factor at the moment of 
horizon crossing proved to be sensitive to the width of the initial pulse. Specifically, by decreasing the width of the 
pulse the amplification factor can drop down by an order of magnitude. This can be understood in terms of the fact 
that parametric amplification strongly depends on the frequency (see discussion in the next Section). 




FIG. 7. Setting an upper limit for the amplification factor: We plot the percentage fractional amplification as a function of 
time. Three difi'erent initial radii are shown: _Ri = lOOM (top curve), _Ri = lOOOM (middle curve) and _Ri = 5000M (bottom 
curve). Different stages of the collapse can be distinguished in the following way: The solid curves terminate at the points that 
correspond to horizon crossing by the stellar surface {R = 2M). The dashed curves continue down to i? ~ 0.25M. The stars 
denote the times where the star has a radius R = 50M. 



In Fig. H we show snapshots of the field itself inside the "mirror" star. The field evolves in a fashion similar to a 
standing wave (there are two main pulses repeatedly bouncing at the left and right boundaries). It is only at the very 
late stages (77 w tt) that amplification can be clearly seen, in agreement with the previous energy considerations. 
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FIG. 8. Diflterent "snapshots" of scalar field propagating inside a "mirror" star. The initial radius of the star is Ri = lOOA/. 
The initial pulse is plotted with a solid curve, the field at 77 = O.Ttt (at which moment has reached the maximum value) is 
denoted by crosses, and finally, the strongly amplified field at 77 = 0.9987r is denoted by the dashed curve. 



V. A PHYSICAL INSIGHT 



In this Section we give some physical arguments which will clarify the apparent inefficiency of parametric amplifi- 
cation. Let us discuss first what happens in the case of cosmological perturbations. There, we have two parameters 
which regulate the strength of parametric amplification: the wavelength A of the perturbation and the Hubble radius 
i?H = a/d. It is well known [ pT[ that amplification is more efficient for large, "superhubble" modes, i.e. A > i?H- This 
is easily understood, as it is a general property of wave propagation in a gravitational field that long wavelengths are 
the ones that are dominantly affected (an example is provided by Fig |6|). 

In the case of a collapsing "mirror" star, there is a third lengthscale entering the problem, namely, the size Xo of 
the star. As evident from the solution ( |l9|) , the field inside the star is made of the discrete spectrum An = 2xo/n. 
Hence, the maximum wavelength is just equal to the diameter of the star. For amplification to be significant, we 
should demand that, 



1 _ 1 
Xo > ~ 2tan(7;/2) ' 

Moreover, there is another restriction set by the time required for the star to cross its Schwarzschild radius, 
condition R > 2M, with the help of eqns. (|4|),(0), gives. 



Xo < 



TT — 77 



(23) 



The 



(24) 



Conditions (|23| ) and (^) are plotted together in Fig. ^ for various values of Xo, or equivalently, of Ri/M. For a given 
choice of Ri/M there is an amplification "window" where there is at least one superhubble mode present. At this 
point we also need some input from the energy expression ( P2| ) , which says that the time-derivative of the potential is 
a crucial factor for amplification. Since i9^F ~ sin(77/2)/ 008^(77/2) we could argue that amplification should become 
stronger as 77 — > tt. In the same limit, and according to Fig. 0, the amplification window shrinks considerably. In 
effect these two counter-balancing factors produce an almost uniform total amplification factor as it was shown in 
Fig. 1^. On the other hand, for the case of pure Friedmann spacetime the corresponding window is significantly larger, 
which explains why amplification can be much stronger (as illustrated in Fig. |l]). 



11 



1.5 



0.5 



-1.5 



2 -2.0 

















r,=io'm 






Ri=10'l\/I 




\ 

X ^ 
\ \ 


r.=io''m 




\l 
\l 
\l 

U 



-2.5 - 



-3.0 



-3.5 



0.0 



0.2 



0.4 



0.6 



0.8 



1.0 



FIG. 9. The amplification "window". The thick sohd hne corresponds to one half the "Hubble radius" Rh as function of 
time, and each thin solid line represents the comoving radius of the star Xo{Ri)- The thin dashed line is plotted for Xo = ti" 
which corresponds to a pure closed Friedmann spacetime. Finally, the thick dashed line is defined by eqn. (M) and represents 
the moment of time that the Schwarzschild horizon is crossed, for the different R\/M. 



VI. CONCLUDING DISCUSSION 



We have presented the results of an investigation on the possibility of having parametrically amplified perturbations 
during gravitational collapse. We have argued, based on the similarity between the time-varying gravitational field 
in the interior of a spherically collapsing star and the field of the closed Friedmann cosmological model, why this 
process will be present, at least in principle. Indeed, our investigation revealed that the process is present, and in 
order to obtain a quantitive estimate of the effect, we have constructed a simplified collapse model (which captures 
the physics relevant to the problem) and studied time-evolutions of a scalar field placed inside the star. We have 
found that, unless the star has shrunk to a radius smaller than the Schwarzschild limit, the scalar field escapes from 
the star in a very short timescale for amplification to become important. We have also quantified the amplification in 
terms of the field energy, and set an upper limit (for scalar perturbations) by considering a perfectly reflecting star, 
i.e. a star that keeps the field trapped in its interior for the entire collapse event (this is what we called a "mirror" 
star). Unfortunately, we found that even under these very favourable conditions, amplification is only of order ^ 1% 
when the stellar surface is about to cross the star's Schwarzschild radius. It is only when the radius has shrunk to 
a much lower value that amplification becomes notable (for example, there is a 10% amplification in the energy 
when R = 0.25M). Moreover, these results remain almost invariant for a wide range of initial radii RJM . We have 
given a simple physical argument, based on the comparison between the stellar radius, the Hubble radius a/d and the 
moment at which the star crosses its Schwarzschild radius, to "explain" the low efficiency of the process. 

The results presented in this paper should encourage one to adopt a pessimistic point of view regarding the astro- 
physical importance of parametric amplification during gravitational collapse. However, this statement is far from 
being definite and must be viewed with some caution. We have considered one of the simplest possible models to 
describe a collapsing star. One could naturally ask how the above conclusion may change when more realistic models 
are adopted (which would incorporate features like pressure gradient, rotation etc.). In all these models, there is 
still present a time- varying background gravitational field which could pump energy into any kind of perturbations 
inside the star. There is no way at this point to say whether the amplification mechanism would operate with the 
same efficiency as in pressureless spherical collapse. Some steps towards including rotational effects could be taken by 
considering a slowly rotating 0-S collapse model |17|. There are also other factors that should be taken into account 
in realistic collapse. For example, it is well established that an inhomogeneous collapsing star can easily give 
birth to a naked singularity instead of a black hole. Under such conditions, our results suggest that strongly amplified 
fields could be generated while the stellar radius is < M , and subsequently escape to infinity. In our view, and despite 
the "negative" results of this study, further (and more detailed) work is needed to give a more complete answer on 
the astrophysical significance of parametric amplification. 
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APPENDIX A: ANALYTIC SOLUTION FOR THE "MIRROR" STAR 



In this Appendix we derive the analytic solution used to describe a scalar field inside the "mirror" star of Section 
IV. We rewrite eqn. (|l^), 

and perform separation of variables: 

We then get the following pair of equations for H and X, 



1 



_ cos-^(r7/2) 



sin^(x) 



H = , 
X ^0 . 



where v is the separation constant. The solution of these equations is, 

X^X) = ciu,l)^^^x)P;!r.hcos{x)], 

2 

H^{ri) ^ ci{v)ui{v,ri) +02(1^)^2(2^,77), 
where c, ci , C2 are constants (with respect to 77 and x) 1 PH [A ^'^'^ associated Legendre functions and 



(Al) 

(A2) 
(A3) 

(A4) 
(A5) 



ui — cos[j/(7r — "n)] — — tan(r7/2) sin[i/(7r — 77)], 
U2 = 2vsnY[v{-K — 77)] + tan(r//2) cos[i^(7r — rj)]. 



(A6) 
(A7) 



The above solutions already incorporate the "left" boundary condition ip{ri,0) = 0. 

At this point, we restrict our attention to the monopole case {£ = 0) only. The time dependence of the field ip (as 
given by eqn. ( [A2| ) is the same irrespective of £, so consideration of ^ ^ will just overcomplicate the analysis without 
providing any new physical information. 

We next impose the "right" boundary condition, 1^(1], Xo) — 0. For £ = we find the following eigenvalues and 
eigenfunctions: 



f„ = — n, for n = ±1, ±2, . . . 

Xo 



V'n = [Cl(l^„)"l(^^n,77) +C2(7^„)u2(7^„,7;)]sin X 

The full field ip can be presented as an infinite sum of these eigenfunctions. The two remaining constants ci,C2 can 
be fixed by the use of initial conditions (assumed static here). The final expression for ip is. 



V'(??,x) = — VC'„ 



nn 
cos [ — rj 

Xo 



Xo . fV\ ■ T^n 

tan — sm — 77 

27rr7. V2/ Vxo 



■nn 
sm [ — X 

Xo . 



where 



Cn= I dx ipoix) sin ( —X ) ■ 
/o VXo 



(A8) 



(A9) 



This solution exhibits a standing wave-like behaviour. Moreover, by setting Xo = it can be used to describe scalar 
wave propagation in a closed Friedmann spacetime. When applying (A8) in practice, we found that the n-sum is 
rapidly convergent (typically, no more than 20-25 terms are required). 
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APPENDIX B: ENERGY "CONSERVATION" LAW FOR THE KLEIN-GORDON EQUATION 

The "master" equation that governs the field -0 is a Klein-Gordon equation. Here we intend to derive a simple 
"energy conservation" expression for this type of equations when homogeneous boundary conditions are assumed for 
the field at the right and left endpoints. 

Consider again eqn. ( p^ ) 

a2^-a^^ + y(?7,x)^-o, (bi) 

where V{ri, x) is an effective potential which can generally have both spatial and time dependence. Suppose that the 
field is confined in the reg ion < x < Xo and ip{ri,0) = '4'iViXo) = (as in the model of the "mirror" star of Section 
IV). Multiplying eqn. (Bl) with 29^^ and making a trivial rearrangement we obtain 

This can also be written as 

a„ [{dri^jf + [d^i^f + v^?] - ij^d,jV - 2d^{dr,i^d^^) = o. (B2) 



Integrating eqn. (B2) from to xo and taking into account the boundary conditions we get 

rxo rxo 

dn / [{dv^f + {dx^f + V^^'] dx = / dxVd^V. (B3) 
Jo Jo 

We define the quantity 

1 rxo 

as the field's total energy (this can be compared with the energy expression for a vibrating string). We have then 
shown that, 

E{t^) - E(Q) r r dri'dx (B5) 
^ Jo Jo 

In the "usual" case of time-independent potentials, this equation just says that the field energy is conserved. This 
property is spoiled when the potential is allowed to vary in time due to the coupling of the field and the time-derivative 
of the potential. This simple relation is an alternative way to view amplification of fields by time-dependent potentials. 
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